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ABSTRACT 

The results of a coordinated space-based photometric and ground-based spectro- 
scopic observing campaign on the enigmatic 7-ray binary LS 5039 are reported. Sixteen 
days of observations from the MOST satellite have been combined with high- resolution 
optical echelle spectroscopy from the 2.3m ANU Telescope in Siding Spring, Australia. 
These observations were used to measure the orbital parameters of the binary and to 
study the properties of stellar wind from the O primary. We found that any broad- 
band optical photometric variability at the orbital period is below the 2 mmag level, 
supporting the scenario that the orbital eccentricity of the system is near the 
0.24±0.08 value implied by our spectroscopy, which is lower than values 
previously obtained by other workers. The low amplitude optical variabil- 
ity also implies the component masses are at the higher end of estimates 
based on the primary's 06.5V((f)) spectral type with a primary mass of 
^26 M and a mass for the compact star of at least 1.8 M . The mass loss 
rate from the O primary was determined to be 3.7 to 4.8 xl0~ 7 M Q yr -1 . 
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* Based on data from the MOST satellite, a Canadian Space 
Agency mission, jointly operated by Microsat Systems Canada 
Inc. (MSCI, formerly the space division of Dynacon Inc.), the 
University of Toronto Institute for Aerospace Studies and the 
University of British Columbia, with the assistance of the Uni- 
versity of Vienna. 



1 INTRODUCTION 

LS 5039 (V479 Set) is the optical cou nterpart of the pe - 
culiar X-ray source RX J1826.2-1450 (|Motch et al.lll997l ). 
The system has been observed in the radio, optical/IR, UV, 
X-ray and 7-ray wavelengths. A ra dio counterpart ha s been 
identified in VLA observations by iMarti' et all (|l998l ). The 
emission appeared to be persistent and was non-thermal. 
A pair of symmetric radio features were later found as- 
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sociated with the source (core) in VLBA observations by 
iParedes et all ||200(ll2002h . which were interpreted as emis- 
sion from two opposite relativistic jets. The system was also 
identified with a very high energy (VHE) 7-ray source found 
in the CGRO/EGKET (Paredes et al. 2000) and HESS 
(|Aharonian et all [2005a) surveys. The detection of MeV- 
GeV/TeV emission placed LS 5039/RX J1826. 2-1450 into 
a class of unusual high-energy objects, the 7-ray binaries. 
So far, only a handful of 7-ray binaries are know n. The 
others ar e PSR B1259-63 llAharonian et all l2005bl ). LS I 
+ 61 303 dAlbert et allkooel ; I Acciari et alj|2008l). Cygn us X- 
1 (| Albert et al.ll2007l ). Cygnus X-3 (iTavani et al.ll2009l). an d 
the recent candidate HESS J0632+057 (|fflnton et al.ll2009l) . 

LS 5039/RX J1826. 2-1450 (hereafter, referred 
to as LS 5039) has been classified as a high- 
mass X-ray binary in the ca tal ogue compiled by 
iLiu. van Paradiis fc van den Heuvell J2006l). The dis- 
tance to the source is w 2.5kpc (|Casares et al.l l2005h . 
Its primary is a bright (V = 11.2 ) massive O star 
(|Clark et al.ll200ll ; [McSwain et alj|200ll ), and its secondary 
is a compact star. It is still a matter of debate whether 
the compact star in LS 5039 is a black hole or a neu- 
tron star. UV and optical spectro scopy established tha t 
the primary is a 06.5V((f)) star (|McSwain et all l2004h . 
The presence of P Cygni profiles in the UV N V A 1240 
and C IV A1550 lines (|McSwain et all \2004 ) indicates 
a strong wind outflow, whose rate has been estimated 
to b e ~ 10 -7 M0yr _1 or even higher (|McSwain et al.l 
20041). From the radial velocities of the Ha and He lines, 



McSwain et al.l (|2004h (hereafter M04) obtained an orbital 



period P = (4.4267 ± 0.0005)d and derived the mass 
function f{m) = (0.0017 ± O.OOO5)M , the orbital ec- 
centricity e = 0.48 ± 0.06 and aisini = (1.36 ± 0.13)R Q , 
where a± is the semi-major axis of the primary's orbit 
and i is orbital inclination. They argued that the mass of 
the primary O star is in the range (20 — 35)Mq, and the 
compact s econdary is a neutr on star which has a mass 
« 1.4M . ICasares et all (|2005l ) (hereafter C05) conducted 
a comprehensive analysis of the optical H Balmer and Hel 
and Hell lines and obtained P = (3.90603 ± 0.00017)d, 
e = 0.35 ± 0.04, f(m) = (0.0053 ± 0.0009)M Q and 
aisini = (1.42 ± O.O7)R . They derived that the orbital 
inclination i = (24.9 ± 2.8)°, the primary O star has 
a mass Mi = 22.91J '9^0 an d the compact star has a 
mass M2 — 3.7^' qMq. The orbital period obtained by 
C05 is shorter than that obtained by M04. It is however 
consistent with the modu l ations observed in th e X-rays 
|Bosch-Ramon et al.l 120051; iTakahashi et al.l 120091) and 7- 



rays at GeV ( Abdo et al.boOS ) a nd TeV 1 Aharonian et all 
120061 ) energies. lAragona et all (|2009l ) (hereafter A09) 
revisited radial velocity (RV) measurements from optical 
spectra, confirming the 3.9d orbital period. Their re- 
fined orbital parameters are P = (3.90608 ± 0.00010)d, 
e = 0.337 ± 0.036, f(m) = (0.00261 ± O.OOO36)M and 
ai sin i = (1.435 ± 0.066)R Q . 

One major question about LS 5039 is the nature 
of the compact object. The orbital parameters de- 
rived by C05 clearly indicate that the compact 
star has a mass exceeding 3.0 Mg, the usually 
accepted upper limit of m asses of neutron stars 
(|Lattimer L Prakashl |2007| ). implying that LS 5039 
is a candidate black-hole high-mass X-ray binary. 



Black-hole binaries with a massive O-type donor star are 
very rare, partly because of the extremely short life-spans 
of such systems. To date, Cyg X-l is the only known black- 
hole X-ray binary in the Milky Way with a massive O donor 
star. The presence of a black hole or not in LS 5039 would 
therefore have significant implications not only on the for- 
mation of black-hole high-mass X-ray binaries and the pop- 
ulation of such systems in the Milky Way, but also on how 
very high-en ergy (TeV) emission is produced in 7-ray bina- 
ries (see e.g. lAraudo. Bosch-Ramon fc Romero! (|2009l )). 

The issue of the nature of the compact star asso- 
ciated with LS 5039 is far from being settled. While 
the spectroscopic observations of C05 favor the 
black-hole scenario, some workers have argued for 
a neutron-sta r scenario involving a non-accreting 
young pulsar dMartocchia et al.l [20051 : iDubusj 120061: 



Sierpowska-Bartosik fc Torresl 2007|: Cerutti et al.l 



2008h . Observationallv. iRibo et al.l (|200a )~ found that 
changes in the mas morphology of the radio images were 
difficult to reconcile with the micro-quasar (black-hole bi- 
nary) scenario in which the radio emission originates from 
expanding plasmons. The situation is further complicated 
by the fact that t he temporal and spectral behaviour of 
the X - ray emission ( Bosch- Ramon et al. 2007; Bosch- Ramon! 
120091 ; ITakahashi et all 120091 ; iKishishita et all |2009[) and 
properties of the TeV e mission jBosch-Ramon et al.l 120081 ; 
iKhangulvan et al.ll2008l ; I Abdo et al.ll2009l ) indicate that the 
high-energy radiation might originate from regions far out- 
side the binary orbit of LS 5039. 

We have been observing LS 5039 photometrically with 
ground-based telescopes at optical and near IR wavelengths 
to search for orbital modulation, but orbital variations have 
not been detected in those data (to be reported elsewhere). 
However, we have detected variations at the 20 mmag level, 
especially in the Jc-band. This prompted us to take a closer 
look at the optical photometric variations of LS 5039 from 
space with the Canadian Microvariability and Oscillations 
of Stars {MOST ) satellite. Our MOST observations were 
made in July 2009, simultaneously with ground-based opti- 
cal spectroscopic observations from the Australian National 
University (ANU) 2.3m Telescope at Siding Spring Obser- 
vatory, Australia. Here we report on the main findings, es- 
pecially those concerning the masses of the component stars 
in LS 5039 and the nature of the compact star. 

This paper is organized as follows. In Section [2] we de- 
scribe the photometric observations with the MOST space 
telescope and the spectroscopic observations with the ANU 
2.3m Telescope, which were made simultaneously. In Sec- 
tion [3] we describe the data analysis, present the results on 
the system's orbital and other parameters, and discuss the 
implications of the nature of the compact object and the 
properties of the wind/outflow from the primary star. The 
conclusions are summarized in Section 2] 



2 SIMULTANEOUS PHOTOMETRIC AND 
SPECTROSCOPIC OBSERVATIONS 

2.1 MOST photometry 

The MOST microsatellite (|Matthews et all Il999l ; 

IWalker et all 120031 ) houses a 15-cm Rumak-Maksutov 
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Figure 1. Continuum normalized spectra of LS 5039, classified 
as type 06.5V((f)), and the 06-7V((f)) spectral template star 
HD 168075. 



telescope feeding a CCD photometer through a custom 
optical broadband (350 - 750 nm) filter. From its polar 
Sun-synchronous orbit (altitude = 820 km; orbital period 
= 101 min), it has a Continuous Viewing Zone (CVZ) 
about 54° wide within which it can monitor target fields 
for up to 2 months without interruption. Targets brighter 
than V ~ 6 are observed in Fabry Imaging mode; fainter 
targets (like LS 5039) are observed in Direct Imaging mode, 
similar to standard CCD photometry with a ground-based 
instrument. The photometry is non-differential, but given 
the orbit, thermal and design characteristics of MOST, 



experience has shown that it is a very photometrically stable 
platform even over long timescales (with repeatability of 
the mean instrumental flux from a non-variable target of 
the brightness of LS 5039 to within about 1 mmag). 

The science exposures on the MOST Science CCD take 
place at the same time as the guide star exposures for satel- 
lite pointing. In the case of the LS 5039 observations, the 
guide star exposure time (and hence the science target ex- 
posure time) was 3.03 s. To build up Signal-to-Noise Ra- 
tio (SNR), the science exposures were co-added on board 
the satellite in 'stacks' of 10 exposures. Stacks (each 30.3 
s of total integration) were downloaded from the satellite 
consecutively (with no dead time between stacks), giving a 
sampling rate of about twice per minute. 

The LS 5039 field lies outside the MOST CVZ so it 
could not be monitored continuously. This field alternated 
with another MOST Primary Science Target field (the Wolf- 
Rayet star WR 113) during each MOST satellite orbit. Typ- 
ically, we monitored LS 5039 for about 70% of every second 
MOST orbit (or about 70 of every second interval of 101 
min) between 2009 July 7 - 23. 

Due to the observing season and the location of the 
LS 5039 field relative to the illuminated limb of the Earth, 
scattered earthshine was high during the beginning and end 
of each 'visit' to the LS 5039 field. We truncated about 15 
- 20 min of each orbit from the original data to preserve 
the highest photometric precision. We also filtered outliers 
caused by cosmic ray hits, especially during passages of the 
satellite through the South Atlantic Anomaly (SAA). How- 
ever, the resulting light curve has an effective duty cycle of 
about 50% (with gaps at intervals of about 101 min) and it 
still sampled the important timescales in the LS 5039 system 
thoroughly. 



2.2 Echelle spectroscopy 

The spectroscopic observations were carried out on four 
nights between 2009 July 8-11 (during the MOST observa- 
tions) and on three nights between August 1-3, using the 
ANU 2.3m Telescope with an echelle spectrograph. In total, 
118 spectra were obtained that cover almost 40 hours with 
nearly uniform sampling of the whole orbit. The integration 
times were between 900 - 1200 s and the spectra covered the 
whole visual range between 3900 A and 6720 A. The nomi- 
nal spectral resolving power was A/AA ~ 23,000 at the Ha 
line, with a typical SNR per extracted pixel of about 100 for 
combined spectra that represented one-hour integrations. 

All data were reduced with standard IRAffj tasks, in- 
cluding bias and flat-field corrections, cosmic ray removal, 
extraction of the 27 individual orders of the echelle spectra, 
wavelength calibration, and continuum normalization. ThAr 
spectral lamp exposures were regularly taken before and af- 
ter every object spectrum to monitor the wavelength shifts 
of the CCD spectra. We also obtained spectra for the tel- 
luric standard HP 177724 and the Q6 -7V((f)) spectral tem- 
plate HD 168075 (jDufton et alj|2006l ). Typical continuum- 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under co-operative agreement with 
the National Science Foundation. 
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Figure 2. Top: Radial velocities based on H Balmer and He II 
lines. Bottom: The best-fitting curve to radial velocities of He II 
lines. 

normalized spectra for LS 5039 and HD 168075 are shown 
in Fig. [T] 



3 DATA ANALYSIS AND RESULTS 
3.1 Radial velocities 

To measure radial velocities we first generated one hour long 
averaged spectra to get higher SNR - one hour corresponds 
to 0.01 orbital phase, hence negligible phase smearing ap- 
pears in the phased radial velocity data. Radial velocities 
of the H I, He I and He II lines were determined by fitting 
two-component functions (a concentric sum of Gaussian and 
Lorentzian functions) and calculating the shift of the cen- 
troids to laboratory wavelengths. To check our method we 
also determined the velocities of some interstellar lines (Ca 
II K, Na I D, DIB A5780, DIB A6613) for each night. We ap- 
plied heliocentric corrections to the radial velocity for each 
line at each time. The data were phased with an orbital 
period of 3.906 d (C05). The success in generating a RV 
curve with this phase folding supports the orbital period of 
LS 5039 as being that obtained by C05 instead of the longer 
4.4 d period obtained by M04. 

To generate the final RV diagrams we used averaged 
velocities of the H I (Ha, H/3, H 7 , H<5, A3835), He I (A4471, 



A5875) and He II (A4200, A4686, A5411) lines; there are sev- 
eral other H and He lines in the wavelength region of our 
spectra (e.g. the frequently used He II A4542 line), but they 
were too noisy or blended to use for velocity determination. 
The resulting RV diagrams are shown in Fig. [2] top. The 
error bars shown in Fig. [5] (typically ±10 to ±15 km s" 1 
in magnitude) represent plus and minus one standard devi- 
ation of the measured velocities of two to five lines at dif- 
ferent wavelengths. The relatively large errors were caused 
partly by observational noise and partly by the high rota- 
tional velocity of the O star (v Iot smi — 113 ± 8 km s _1 , 
C05). 

We detected systematic increasing blueshifts from the 
He II lines to the He I lines to the H Balmer lines. C05 have 
also described this effect, but found its degree to be smaller 
(they obtained an ~8 km s _1 shift between the average RVs 
of the He II and H Balmer lines, as opposed to our 20 km 
s -1 value). We believe that the RV shift is due to the con- 



tamin ation of wind for the He I and H I lines (see lPuls et alj 
(1996)). and so we did not use those two atomic line sets to 
constrain the orbit. 

Additionally, redshifted satellite absorptions were found 
in the Ca II K and Na I Dl and D2 interstellar absorp- 
tion lines with a radial velocity around ±60 km s _1 (more 
precisely, ±58.4 ± 2.2 km s _1 by Ca K, ±62.9 ± 2.3 km 
s" 1 by Na Dl and ±61.9 ± 1.7 km s" 1 by Na D2). These 
satellite lines may belong to a formerly unknown Galac- 
tic Intermediate Velocity Clo ud [IVC; see the review by 
IWakker fc van Woerder] (| 19971 )]. 



3.2 Orbital and system parameters 

In the following analysis we adopted T c s = 39,000 ± 1000 
K and log g = 3.85 ± 0.10 for the O primary (C05). We 
measured the equivalent widths (EW) of several interstel- 
lar lines to estimate the interstellar reddening (see Section 
13.41 for EW measurements of other lines). Using Na I Dl 
[E W = 0.70 ± 0.02 A with the relation to reddening as 
per iMunari fc Z witter! i|l997h ], DIB A5780 and DIB A6613 
[EW = 0.55 ± 0.05 A and 0.18 ± 0.02 A, r espect ively with 
the relation to reddening as per lCox et al! (|2005l )] lines, we 
found E(B — V) = 1.2 ± 0. 1, which is i n agr eement with 
previous results [1.25 to 1.35. lRib6 et al! (|2002f l and 1.28 ± 
0.02, M04]. Based on this agreement, we adopted the values 
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d = 2.5 ± 0.1 kpc, Mo = 22.9 
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velocity points and did not use previous RV measurements 
made by others. Our data set therefore represents the 
highest resolution, homogenous spectral data set yet used 
to obtain an orbital solution for LS 5039. Given the lack 
of X-ray eclipses we could not use the special mode in 
WD code developed for modeling the orbits of X-ray 
binaries. Therefore, we analysed LS 5039 as a single line 
spectroscopic binary without any (X-ray) light curves. This 
limitation allowed us to determine only the mass function, 



ftp: / /ftp. astro. ufl.edu/pub/wilson 
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Table 1. Orbital parameters of LS 5039. 
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Parameter 



C05 (He II) 



A09 



This paper (He II) 



T (HJD-2450000) 1943.09 ± 0.10 2825.99 ± 0.05 5017.08 ± 0.06 



Porb (d) 

e 

W n 

V 1 (km s" 1 ) 
A'i (km s- 1 ) 
a\ sin i (Rq) 
f(m) (M Q ) 



3.90603 
0.35 ± 0.04 
225.8 ± 3.3 
17.2 ± 0.7 
25.2 ± 1.4 
1.82 ± 0.10 



3.90608 
0.34 ± 0.04 
236.0 ± 5.8 
4.0 ± 0.3 
19.7 ± 0.9 
1.44 ± 0.07 



3.906 (adopted) 
0.24 ± 0.08 
237.3 ± 21.8 
3.9 ± 1.3 
23.6 ± 4.0 
1.77 ± 0.15 



0.0053 ± 0.0009 0.0026 ± 0.0004 0.0049 ± 0.0006 



RMS of fit (km s" 1 ) 9.1 



7.1 



6.2 



f(m), as a function of different inclination angles, i. We 
could not determine the exact values of the inclination, i, 
and the mass ratio, q. However, knowledge of the primary's 
mass, Mq, allowed us to narrow the possible parameter 
space. 

As described previously, we adopted the 3.906 d period 
as a fixed parameter during the modelling. The computed 
value of To = HJD 2455017.08 was used as the epoch of peri- 
astron. The computed orbital parameters are given in Table 
[1] where they are compared to the results of C05 and A09. 
The RV curve implied by the solution is shown in Fig. [2] bot- 
tom. Our values and the ones published by C05 are based 
only on velocity points from the He II lines, while A09 ap- 
plied the velocities of every available H, He I, and He II line, 
a process which combines lines from two different sources on 
and near the O star as mentioned in Section \3. II 

In general, our computed orbital parameters are close to 
earlier solutions, but there are some differences. The value of 
the computed systemic radial velocity, V 11 was significantly 
higher (by 15 to 20 km s _1 ) for each line type (H I, He I 
and He II) in C05 than what we found. The possibility of a 
real change in the system RV over a few years is very small, 
so the cause of the difference is likely due to differences in 
data analyses. 

3.3 Mass of the compact object 

One of the main goals of our investigation was to obtain 
stronger constraints on the mass of the non-stellar com- 
panion. C05 executed detailed light curve simulations using 
their orbital solution for LS 5039; they found that if the in- 
clination angle is around 30°, then photometric variability 
caused by the distortion of the primary should be of the or- 
der of 0.01 mag near periastron. If the change in brightness 
is 0.01 mag or less then the inclination is 30° or less. An 
inclination less than 30°, in turn, implies that the mass of 
the compact object is too high (> 3.0 Mq, C05) to be a 
neutron star. Lacking the necessary 2-3 mmag photometric 
accuracy required, they could not check their scenario. 

As described in Section 13.21 our new analysis of the or- 
bital parameters of LS 5039 is based on an independent ho- 
mogeneous radial velocity data set. Our results are in good 
agreement with the one presented by C05 - in particu- 
lar, the values of mass functions, the implications of 
which are shown in Fig. [3j agree within the uncer- 
tainties. 

Photometric data from the MOST satellite (6649 in- 
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Figure 3. Mass constraints of the two components in LS 5039 
for different inclinations from the measured mass function. The 
vertical lines arc the limits of the O-type stellar companion mass 
adopted from C05. The gray region represents the possible values 
of the compact object given that the orbital inclination is less 
than 30° . 

dividual brightness measurements through 16 days between 
July 7 and July 23) did not show any variability greater than 
a few mmag. To quantify the frequency content of the light 
curve, we performed a perio d analysis of the fu ll dataset 
using the Period04 software (|Lenz fc Bregei|[2003 ). The re- 
sulting frequency spectrum does not contain any significant 
peak with an amplitude greater than 0.002 mag. Moreover, 
the orbital period does not jump out of the noise either. 
As a different approach, we phased the light curve with the 
spectroscopic period and epoch and used 0.03 phase bins 
to reveal if there is any optical variability due to orbital 
motion. Fig. [4] shows the phase-binned folded light 
curve, which indicates a possible variability at the 
level of 2 mmag, with an apparent broad minimum 
at phase tp ~ 0.7 - 0.8. Although the SNR of the binned 
light curve is low for detailed speculations, we compared the 
shape of the curve with the EW variations of hydrogen and 
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orbital phase 

Figure 4. Phase- binned folded light curve from the MOST obser- 
vations. The error bars represent 99% confidence intervals of the 
mean values computed from 3<r/y/n where a is the standard 
deviation and n is the number of data points in a given 
phase bin. The MOST data originally represent fractional 
change from the mean in terms of flux so a change of sign 
gives the A mag values shown here that represent the 
magnitude difference from the mean. This conversion of 
the MOST data to A mag allows us to directly compare 
the data to the simulations shown in Fig. QjJ 
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Figure 5. Model light curves. All models have the mass 
function fixed to 0.0049 Mq. Mq is the mass of the pri- 
mary, M x is the mass of the compact companion. The 
amplitude of the light curves decrease with increasing 
system mass and decreasing eccentricity. 



helium lines (see Section I3.4[) and with the results of 
our own light curve simulations. 

To check and understand the conclusions of the 
C05 light curve simulations, we used the WD code to 
do our own simulations. Two sets of simulated light 
curves were computed: for the first set the mass of 
the compact object was fixed to 3.0 Mq; for the 



second set the mass of the compact object was de- 
termined from the mass function derived from our 
radial velocity data as described in Section I3.2I In 
each set, the mass of the primary was set to either 20 
or 26 Mq, the eccentricity was set to either 0.35 or 
0.24, the V band was used to approximate the MOST 
band and, inclinations of 60°, 30° and 20° were mod- 
elled. The results from the second set of simulations, 
with the mass function fixed to our measured value, 
are shown in Fig. [5j For all simulations the primary 
radi us was set to Rp = 9.5 Rq a nd a linear cosine 
law (|Wilson &i van Hammdl2003h was used for limb 
darkening effects. The compact object was handled 
as a point source. 

The first set of simulations, with the mass of the 
compact object fixed at 3.0 Mq, gave light curves 
that followed the pattern described by C05: the am- 
plitude of the light curve decreased as the inclination 
decreased. The results with the fixed mass function 
simulations were different, as seen in Fig. Q3 the am- 
plitude of the light curve decreased with increasing 
total system mass and decreasing eccentricity but 
did not decrease with decreasing inclination. For 
fixed system mass and eccentricity the inclination 
diagnostic is better given by the light curve shape, 
especially the dip near phase 0, and not by its am- 
plitude. A formal curve fit of the MOST data to the 
light curves of the Mo = 26 Mq, e = 0.24 case gave the 
best fit with an inclination of 60°. However this re- 
sult must be viewed as preliminary, at best, because 
the amplitude difference between the different incli- 
nation light curves is very small, the binned MOST 
data have considerable scatter at the millimag level, 
and the effect of the stellar wind and accretion on 
the light curve are not modelled in the WD code. 

We note that some studies suggest the mass of an Q6.5V 
type star to be around 28 - 29 M Q (jMartins et all I2005T ) 
rather than 20 — 26 Mq used in C05. In the case of 
a heavier primary, the mass of the compact object 
would also be larger and the M x masses shown in 
Fig. [5] would be greater, strengthening the argument 
for the black hole nature of the secondary compact 
star, even at relatively higher inclinations. 

3.4 Stellar wind 

To infer the mass loss rate of the O-type star and the prop- 
erties of the circumstellar matter, we determined the equiv- 
alent widths of each H and He line and also their variability 
during the orbit, which could be good indicators of physical 
processes taking place in the stellar wind. As mentioned in 
Section 3.2, we also measured the EWs of some other lines 
to check the value of interstellar reddening given in the liter- 
ature. We made the measurements on the average of spectra 
summed to one hour long exposure times, but we plotted the 
daily average values of the EWs on final diagrams, shown in 
Fig. [HI to better see the trends. 

For the Ha line we found that EW changes from 2.50 
to 2.85 A over an orbital phase (see Fig. The average 
value of 2.70 A agrees very well, within the uncertainties, 
with the result of C05 (2.8 ± 0.1 A), except they found 
that the value was stable during the orbit of the binary 
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for the O-type star (C05); a terminal wind velocity of Voo 
= 2440 ± 190 km s _1 with a wind velocity law exponent 
of P = 0.8 (M04). We found that the mass- loss rate of the 
stellar companion is around 3.7 x 10~ 7 Mq yr _1 from the 
strongest absorption, which corresponds to the lower limit, 
and 4.8 x 10 -7 Mq yr _1 for the upper limit. These values 
are consistent with the mass loss rates of ~ 3.7 x 10~ 7 Mq 
yr^ 1 for the low state (strong absorption) and ~ 7.5 x 10~ 7 
Mq yr for the high state obtained by C05 (see also M04). 

We found two lines (H/3 and He I A5875) showing signif- 
icant changes during the orbit. The lowest absorption for the 
H/3 line occurs around y? ~ 0.75, and at ip ~ 0.65 for the He 
I line (Fig. [6]), close to the expected phase ip ~ 0.7 when the 
compact object is between us and the stellar companion (in- 
ferior conjunction). The lower EW at inferior conjunc- 
tion implies an increased emission strength likely 
due to the focusing of the stellar wind toward the 
compact object. Models of the wind flow in the sys- 
tem will need to take that focusing into account, es- 
pecially those that model a pulsar wind/stellar wind 
interaction. We carried out a simple correlation analysis 
between these EW changes and the MOST light curve (re- 
binned to the EW bin size), and found the correlation 
coefficients, r, to be 0.52, 0.70 and 0.52 for the Ha, H/3 and 
He I A5875 lines, respectively. This result suggests the possi- 
bility of real modulations over the orbital period, but further 
studies are necessary. 

Note that our data only show smooth orbital modula- 
tion in the Ha, H/3 and He I absorption lines and in the 
emission components of the H Balmer lines. The lines do 
not show evidence of dense clump-like condensates in the 
stellar wind. 



He I 5875 A 



Figure 6. Variability of equivalent widths of Ho, H/3 and He I 
A5875 absorption lines during the orbit. The mean values and the 
given errorbars of points were calculated from results we got from 
three different methods (direct measurement with IRAF splot 
task, and fitting Gaussian or Lorentzian curve to the line profile) 
used to determine EWs of lines. 



(which could possibly be explained by the relatively low res- 
olution of their spectra). Our result is also consistent with 
the EW values measured by others over the las t ten years 
Bosch-Ramon et al.ll2007f ). Using the method of I Puis et al.l 



19961 ) we estimated the mass loss rate from the EW of the 



Ha line. To do these calculations we adopted the following 
parameters: R Q = 9.3±j};£, R© and T cff = 39,000 ± 1000 K 



3.5 Some remarks 

The two main findings from our work are that the orbit has a 
significantly lower eccentricity than previously thought and 
that the total mass of the system is toward the higher 
end of previous estimates, when the variability of 
the MOST data is considered. We were not able to 
resolve the nature of the compact object, between 
neutron star and black hole, in spite of simulations 
by C05 that postulated that a low amplitude light 
curve would be associated with a low inclination and 
a black hole. The eccentricity e « 0.4 obtained by M04 put 
LS 5039 as the system with the most eccentric orbit among 
X-ray binaries with an O donor star. Although the value 
was later shown to be lower (e ~ 0.35, C05 and A09), the 
eccentricity would still be high enough to have some con- 
sequences on the dynamical and radiative properties, and 
on the evolutionary history, of the system. The observed 
TeV spectrum is consistent with inverse Compton 
scattering from electrons distant from the compact 
object either in a pulsar wind/stellar wind shock if 
the compact object is a neutron star, or in a jet if 
the compact object is a black hole. The lack of ac- 
cretion features in the X-ray spectrum favors the 
pulsar model unless the bulk of the accretion power 
is released as kinetic energy in jet o utflow instead of 
therm ally from an accretion disk (|Takahashi et ail 
l2009l l . In either case the keV X-rays are hypothe- 
sised to originate as synchrotron emission from the 
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electrons responsible for the inverse Compton TeV 
emission. This is in contrast to the canonical view 
of X-ray emission from X-ray binaries. In the canon- 
ical view, the keV X-rays from X-ray binaries with an 
star are powered by the accretion of gas either in a wind 
focused by the Roche-lobe of the donor star or captured di- 
rectly from the wind out-flowing from the donor star. In the 
focused wind scenario, an accretion disk is expected to form 
around the compact star, similar to that in systems with 
Roche-lobe filling mass-transfer, as the accreting gas carries 
substantial specific angular momentum. In the direct wind 
capture scenario, provided that the specific angular momen- 
tum of the gas is sufficiently small, an extensive accretion 
disk might not be formed and the accretion inflow would 
be practically radial and resemble that of a Bondi-Hoyle 
flow. A lower orbital eccentricity would give a lower spe- 
cific angular momentum in the wind material that is swept 
up by the accreting star. So far observati ons [e.g. X-ray ob- 
servations by Chandra and XMM -Newton (|Martocchia et all 
l2005l ; lB~osch- Ramon et alj|2007l )] have not shown evidence of 
an accretion disk in LS 5039. If LS 5039 has an extremely 
high orbital eccentricity as originally measured, the lack of 
an accretion disk around its compact star would need cer- 
tain non-trivial explanations (cf. the situation in the Be X- 
ray binaries). The lower eccentricity may ease the situation 
somewhat, allowing radial gas inflow for a substantial dis- 
tance before reaching the compact accretor. Then, if there 
is an accretion disk, it would not be expected to be a large 
and dense disk because such a disk would generate signif- 
icant thermal soft keV X-rays, which are not detected in 
the Chandra, XMM-Newton and other X-ray observations. 
Our smaller value of eccentricity reduces the discrepancy be- 
tween the low observed X-ray variability and what would be 
expect ed from wind ac c retion with Bondi-Hoyle lik e radia l 
inflow [Bondi fc Hovlel { 19441 ). iBosch-Ramon et all (|2005h . 
C05]. 

Known black-hole high-mass X-ray binaries tend to 
have small orbital eccentricity, unlike the X-ray pulsar s (see 
(|Liu. van Paradiis fc van den Heuve]|2005l . l200al2007M . X- 
ray binaries with a black hole and an O donor star are very 
rare, and Cyg X-l is the currently only known system in 
the Milky Way. The confirmation of LS 5039 as a black-hole 
high-mass X-ray binary in which a massive O donor star and 
a black hole revolve around each other in an eccentric orbit 
would have important implications on how such systems 
are formed and how massive binaries evolve. In LS 5039, 
the compact star's progenitor, at a certain stage, should be 
more massive than the current O donor star, otherwise it 
would not have evolved to form the compact star. O stars 
have very short life spans (~ a few Myr), thus LS 5039 as an 
X-ray binary must be younger than a few million years. 
This is supported by the fact that LS 5039 has a highly ec- 
centric binary orbit, which has yet to be circularized from 
a presumed recent supernova event. However, it is unclear 
whether the progenitor of the O star or the progenitor of 
the compact star had the larger initial mass, as a rapid 
evolution of the progenitor of the compact star could well 
be triggered by a mass transfer process, in which the system 
was compact enough to allow the progenitor of the current 
O star to overfill its Roche-lobe and transfer material to the 
progenitor of the compact star. ( For more on the evolutio n 
of very massive binaries, see e.g. iDalton fc Sarazi n (1995); 



IVan Bever fc Vanbeverenl (|2003f ); lDionne fc Robert] (|2006l V) 
If LS 5039 was indeed formed through this mass-transfer 
channel, a finding that it contains a black hole among 
the fewer than 200 known high-mass X-ray binaries in the 
Milky Way would imply that a black-hole high-mass X-ray 
might be formed in close massive binaries more easily than 
previously thought. Young stellar clusters in star-forming 
galaxies may well be populated by LS 5039-type sources and 
are potentially Gev-TeV 7-ray sources. 



4 CONCLUSIONS 

Our simultaneous optical photometry from the MOST space 
telescope and high resolution echelle optical spectroscopy 
from the ANU 2.3m Telescope have put further constraints 
on the orbital parameters of the LS 5039 system. In particu- 
lar we obtained a mass function f(m) « 0.0049 ±0.0006 Mg 
and an orbital eccentricity e = 0.24 ± 0.08. The maximum 
photometric variation of LS 5039 in the MOST light curve 
was 2 mmag. The low photometric variation is consis- 
tent with the lower orbital eccentricity of 0.24, as 
opposed to the value of 0.35 found by others, and 
it is consistent with a higher overall system mass. 
However, we cannot determine the inclination on 
the basis of the low photometric amplitude. Our 
light curve simulations imply that the mass of the 
compact object is at least 1.8 TSAq based on an in- 
clination of 60° or less. Our value for the eccentricity of 
0.24±0.08 is a little smaller than previous determinations. 
The lower eccentricity implies that the wind material that 
is captured and falls into the Roche lobe of the compact 
star has a lower specific angular momentum. Thus it may 
not lead to the formation of a large-scale optically thick ac- 
cretion disk, and the accretion inflow is practically radial 
resembling that of a Bondi-Hoyle flow. This radial inflow 
would explain the lack of an accretion disk signature 
in keV data and should be considered in models of 
mass transfer for either a non-accreting pulsar or a 
black hole. Finally from EW measurements of the Ha line, 
we derived that the mass loss rate from the O-type primary 
through stellar wind is 3.7 to 4.8 x 10 -7 Mq yr _1 , similar to 
values obtained by other workers. Our observations do not 
show evidence of dense clumps in the stellar wind. 
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